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Abstract 

We investigate a torsional pendulum containing a ferrofluid that is forced periodically to undergo 
small-amplitude oscillations. A homogeneous magnetic field is applied perpendicular to the pendu- 
lum axis. We give an analytical formula for the ferrofluid-induced "selfenergy" in the pendulum's 
dynamic response function for monodisperse as well as for polydisperse ferrofluids. 
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I. INTRODUCTION 

Real ferrofluids [lj contain magnetic particles of different sizes Q. This polydispersity 
strongly influences the macroscopic magnetic properties of the ferroffuid. We investigate 
here the effect of polydispersity on the dynamic response of a ferrofluid pendulum. 

A torsional pendulum containing a ferrofluid is forced periodically in a homogeneous 
magnetic field H ex t = H ext e x that is applied perpendicular to the pendulum axis e z (see 
fig. HJ). Such a ferrofluid pendulum is used for measuring the rotational viscosity [3]. The 
cylindrical ferrofluid container is here of sufficiently large length to be approximated as 
an infinite long cylinder. We consider rigid-body rotation of the ferrofluid with the time 
dependent angular velocity fi = ipe z as can be realized with the set-up of jjj. The fields H 
and M inside the cylinder are spatially homogeneous and oscillating in time. 



II. EQUATIONS 

First, the Maxwell equations demand that the fields H and M within the ferrofluid are 
related to each other via 

H + AM = H ex * (2.1) 
with N = 1/2 for the infinitely long cylinder. Then we have the torque balance 

Cp = -ufo - r> - 1 + /(*) (2.2) 

with the eigenfrequency u and the damping rate T of the pendulum without field and the 
total moment of inertia G. The magnetic torque reads 

T — —Ho J dV(M x H) t = -VoV(M x W xt ) z , (2.3) 

and f(t) is the external mechanical forcing. 

Finally, we need an equation describing the magnetization dynamics. Here, we consider 
the polydisperse ferrofluid as a mixture of ideal monodisperse paramagnetic fluids. Then 
the resulting magnetization is given by M = ^Mj, where Mj denotes the magnetization 
of the particles with diameter dj. We assume that each Mj obeys a simple Debye relaxation 
dynamics described by 

d t Mj — Q, x Mj = --[Mj - Mf (H)] (2.4) 



February 2, 2008 



3 



We take the equilibrium magnetization to be given by a Langevin function 

Mf (H) = Xj (H)H = Wj C (^^4^) | ( 2 -5) 

with the saturation magnetization of the material M mat and the magnetization distribution 
Wj(dj). Note that the magnetization equations ()2.4j) for the different particle sizes are 
coupled by the internal field H = H eii — iVM. As relaxation rate we combine Brownian 
and Neel relaxation — = \ + \. The relaxation times depend on the particle size by 

t b = W^jridj + 2s) 3 and r J N = /q" 1 exp ( j^pr j with 77 the viscosity, s the thickness of the 
nonmagnetic particle layer, and K the anisotropy constant. 
Altogether we use the following system of equations: 

<p = Q (2.6) 
Q = _ w 2p _ r fi - ^H ext M y + f(t) (2.7) 

Ml = -nMl --[Ml- Xj (H)(H ext - NM X )] (2.8) 
T j 

Mi = tlMi - -M 3 y - -N Xj (H)M y . (2.9) 

III. LINEAR RESPONSE ANALYSIS 

For the equilibrium situation of the unforced pendulum at rest that we denote in the 
following by an index one has ip = Q = M^° = and = Mi q (H ). Furthermore, 
M = E M 4(#o) with H solving the equation H = H ext - NM {H ). 

External forcing with small |/| leads to small deviations of ip, of Q, and of 5H = H — H = 
— iV(M — M ) = — N5N1/2 from the above described equilibrium state. We expand each 
Xj{H) up to linear order in 5H 

Xi (|H + 6H\) = X jo - x' j0 NSM x + 0(5H) 2 . (3.1) 

Here, Xjo = Xj(Ho) and x'jo * s the derivative of Xjo- Then we get the linearized equations 

ip = VI (3.2) 

Q = -u%(p-r n-Ky + f(t) (3.3) 
1 1 



UJ - ^ ^N( Xj o + x' j0 H )x (3.4) 



Vj = Qx °j ~ —Vj - — N XjoV • (3-5) 
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We have intoduced the abbreviations Xj = 5M X /M , x® = M^' /M , yj = 5M y /M and 
x = Ylj x j> V = TltjUj- The strength of the coupling constant between the mechanical 
degrees of freedom ip, Q and the magnetic ones is k = fioH ext MoV/Q . 
For periodic forcing f(t) = fe~ lut we look for solutions in the form 



fit) 

Xj(t) 



V 



/ \ 

<p 
Ci 



-iujt 



(3.6) 



%•(*) J \Vi J 

Inserting the ansatz ()3.6|) into the linearized equations ()3.2|) - (|3.5|) yields 
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(3.7) 
(3.8) 
(3.9) 

(3.10) 



and 



ip — Gf— [ul -lu 2 - iluTo - cjS] 1 / . (3.11) 

The ferrofluid-induced selfenergy £(u;) in the expression for the dynamical response function 
G(u) of the torsional pendulum is 

„0 



IK 



^ 1 - ium / 4^ 1 



TjXj 



lUJTj 



— IUT. 



(3.12) 



3 / 3 

Its imaginary part changes the damping rate T of the pendulum for k — 0, i.e., in zero 
field. The real part shifts the resonance frequency of the pendulum. In the special case of a 
monodisperse ferrofluid on has 



EM 



1 — iujr + Nxt 



(3.13) 



IV. RESULTS 



We evaluated the linear response function G{uS) = <p(u)/f of the pendulum's angular 
deviation amplitude <p(uj) to the applied forcing amplitude / and the selfenergy £(cu) for 
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some experimental parameters from Q]: co>o/27r = 32.7Hz, T = 0.178Hz, V/Q = 20m/kg. 
The cylinder is filled with the ferrofluid APG 933 of FERROTEC. Therefore, we used in 
equation 1)3.13)1 an experimental r = 0.6ms and the experimental M eq (H) shown in fig. El 
These monodisperse results were compared with the expression 1)3.12)1 for the polydisperse 
case for the typical parameter values M mat = 456kA/m, rj = 0.5Pa ■ s, s = 2nm, K = 
AAkJ/m 3 and /o = 10 9 Hz. The contributions w(dj) that enter into the formulas (J2.5j) for 
the susceptibilities Xj are given by a lognormal distribution j^: 

w(dj) = M sat |^ with g(d 3 ) = -=J exp ( J** {4A) 



Fitting the experimental M eq (H) with a sum of Langevin functions (J2.5)) yields M sat = 
1814JL4/m, d = 7nm and a = 1.47 (see fig. |2J). We used here 30 different particle sizes 
from d± = lnm to d^o = 30nm (see fig. EJ)- 

The calculations show the additional damping rate caused by the interaction between 
ferrofluid and external field. An increasing magnetic field leads to smaller amplitudes; in 
polydisperse ferrofluids the amplitude decreases faster [fig. HI and El (a)]. Furthermore, 
one can see a shift of the peak position to higher frequencies oo maxi which is stronger in 
polydisperse ferrofluids [fig. Eland [5] (b)]. 
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FIG. 2: x: Experimental equilibrium magnetization M eq (H) used as input for the monodisperse 
calculations; full line: fit with lognormal contribution. 
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FIG. 3: lognormal contribution w(di) (di = lnm . . . d^Q = 30nm) used as input for the polydisperse 
calculations. 
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FIG. 4: \G\ near the resonance too; x H ext = OkA/m, squares H ext = 5kA/m, circles H' 
WkA/m; filled symbols: polydisperse. 




FIG. 5: Maximum value max\G\ (a) and peak position u> max (b) as a function of external field 
H ext ; full line monodisperse, dashed line polydisperse. 
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FIG. 6: Re(T,) (a) and Im(L) (b) at u = uj max ; full line monodisperse, dashed line polydisperse. 



